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ABSTRACT 

The X-ray emission by hot gas at the centers of clusters of galaxies is commonly modeled 
assuming the existence of steady-state, inhomogeneous cooling flows. We derive the metallicity 
profiles of the intracluster medium expected from such models. The inflowing gas is chemically 
enriched by type la supernovae and stellar mass loss in the outer parts of the central galaxy, which 
may give rise to a substantial metallicity gradient. The amplitude of the expected metallicity 
enhancement towards the cluster center is proportional to the ratio of the central galaxy luminosity 
to the mass inflow rate. The metallicity of the hotter phases is expected to be higher than that 
of the colder, denser phases. The metallicity profile expected for the Centaurus cluster is in good 
agreement with the metallicity gradient recently inferred from ASCA measurements (Fukazawa et 
al. 1994). However, current data do not rule out alternative models where cooling is balanced by 
some heat source. The metallicity gradient does not need to be present in all clusters, depending 
on the recent merging history of the gas around the central cluster galaxy, and on the ratio of the 
stellar mass in the central galaxy to the gas mass in the cooling flow. 

Subject headings: Cooling Flows — Galaxies: Clusters: General — Galaxies: Clusters: Individual 
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1. INTRODUCTION 

X-ray observations in clusters of galaxies have shown that the central gas often has short cooling 
times (Edge et al. 1992), and that the X-ray spectra show an excess of low-energy photons compared 
to single temperature spectra. This excess may be understood as arising from cold gas components at 
temperatures well below the average X-ray temperature. It is usually inferred from these observations 
that the gas is cooling and that a "cooling flow" forms (see, e.g., Fabian 1994 for a review). 

When interpreted in terms of steady-state cooling flow models, the X-ray surface brightness 
profiles in cooling regions yield a mass inflow rate that increases with radius, requiring an 
inhomogeneous medium with phases at different temperatures at every radius, generating a 
distributed mass drop-out as some of the gas cools below X-ray temperatures (Nulsen 1986; Thomas, 
Fabian, & Nulsen 1987). However, it has not been proved that such a multiphase medium can be 
maintained against buoyancy forces; moreover, the presence of a gas inflow has not been directly 
demonstrated, and the fate of the cooled gas remains unknown. A possible alternative is the presence 
of an energy source that can balance cooling, and for which several alternatives have been proposed 
(Tucker & Rosner 1982; Silk et al. 1986; Miller 1986; Pringle 1989). This could lead to a cycle where 
the heated gas could rise buoyantly from the center, while cooled gas could form clouds which would 
fall and re-evaporate, with no need for a net inflow or mass deposition (see Tabor & Binney 1993). 

Recently, a metallicity gradient was observed in the hot gas of the Centaurus cluster (Fukazawa 
et al. 1994). In this Letter we consider the expected metallicity profile in inhomogeneous cooling 
flow models, as well as in alternative models with no net inflow. Although it might seem that in 
cooling flow models one should not expect metallicity gradients since the gas is flowing in from large 
radius, we show that a substantial metallicity gradient may result due to metal enrichment of the 
inflowing gas by the central galaxy. 

2. COOLING FLOW MODELS 

Inhomogeneous cooling flow models (see Nulsen 1986) assume that at all radii there is a 
multiphase medium with a distribution of temperatures. Relative motions of gas in different phases 
due to buoyancy and heat exchange among phases are prevented by a tangled magnetic field. While 
it has not been shown that these assumptions can be valid in reality, they are the only known way to 
reconcile steady-state cooling flow models without energy injection with the observed X-ray profiles 
and spectra. The phases that cool to low temperatures are "deposited" and no longer contribute to 
the average weight of the multiphase medium, so the mass inflow rate M[r) increases with radius. 

Cooling flows generally occur around central cluster galaxies, which often have extended cD 
halos. Evolved stars in these galaxies should inject metals into the hot gas at a rate ([r) proportional 
to the stellar density (neglecting any radial metallicity gradients for the stars). Since metal injection 
from stars is equally likely to occur at any point, we assume that the metallicity in a phase of density 
p is increased at a rate Z{r,p) = C,{r)/ p. Thus, the hotter phases become more metal rich than the 
cooler (and denser) ones. Metals could be produced by type-la supernovae or lost from evolved stars 
with higher metallicity than the gas. The metal-rich ejected gas will have a large kinetic energy 
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(from the supernova explosion or from stellar orbital motions) which will heat the surrounding gas, 
further increasing the concentration of ejected metals in the hot phases. 

The power injected by all the supernovae required to produce an iron abundance enhancement 
AZ in the cooling flow region is Lsn ~ [AZM lMFe)E, where E ~ 10"erg and Mpe Z 0.3M@ 
are the kinetic energy and iron mass produced by a type la supernova (e.g., Woosley & Weaver 
1986). The total luminosity of the hot gas in the cooling flow region is i ~ (5T/2/x)M, where T is 
the temperature of the gas, and fi = 10~'^'*g is its mean molecular weight. Their ratio, Lsn /L ^ 

0. 1(AZ/O.5Z0)(T/3 keV)"-*^, is small for the Centaurus cluster (for which the reference parameters 
were chosen), and probably even smaller for richer, hotter clusters, so that the contribution of 
supernovae to the total energy budget is negligible. Evolved stars will inject gas at their virial 
velocity, which is comparable to the gas thermal motions. In order for the cooling flow model to 
be meaningful, the mass of gas injected has to be a small fraction of the total gas mass present 
(requiring a stellar metallicity much higher than the maximum metallicity of the gas), so that the 
heat budget is again not substantially altered. Thus, we assume in what follows that only metals 
are injected, but the mass and energy budget of the cooling flow are unaffected by this process. 

If we identify a phase by its "initial" density po at a fiducial radius t-q, its metallicity at radius 

r is 

Z{r,po) = Z^+ r ,'^^"'\^"'. , (1) 
Jr P{r',Po)u{r') 

where p[r, po) is the density at radius r of the phase of initial density po, and u[r) is the inflow 
velocity of the gas, assumed to be the same for all phases. Zoo is the metallicity at a large distance 
from the central galaxy, which we assume was laid down before the cluster gas relaxed to its present 
steady state. In principle, the integral should have an upper cutoff at a radius comparable to the 
cooling radius, rcooi (the radius at which the cooling time equals the age of the cluster), beyond 
which the steady-state assumption does not hold. Since most of the metals are injected at smaller 
radii, we will ignore this upper cutoff. 

If the total enclosed mass, the gas density profile, and the gas cooling function in the range 
of interest for the cooling flow are well approximated by power laws, M[r) oc r-'*'", p oc r~^'', and 
A(T) = AqT", then the cooling flow is self-similar (Nulsen 1986; Waxman & Miralda-Escude 1995), 

1. e., the mass inflow rate at radius r of gas in phases with density larger than p can be written as 
tp^r, p) = M[r)I[p/ p[r)], where p[r) is the average gas density at radius r. The total mass inflow 
rate is also a power of radius, M oc t"*, with 6 = 3 — 2\p — (1 — a)[\M — 1)- The cumulative density 
distribution at any given radius is 
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otherwise; 

with v = 56/{2 - a)/[2\p + 3(Am - 1)]. 

The density evolution of a given phase is given implicitly by 
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where rdepi^jp) is the deposition radius of the phase whose density is p at radius r, and Pmin{f) 
and po,min are the densities of the hottest phase at radius r and at the fiducial radius, respectively. 
Thus, for this case eq. (1) can be rewritten as 



Here, we apply the results of § 2 to the Centaurus cluster (A3526), in order to show that a 
substantial metallicity gradient can in fact be expected in a cooling flow model. We chose this 
cluster because its radial metallicity distribution has been measured from spatially resolved X-ray 
spectroscopy by ROSAT (Allen & Fabian 1994) and ASCA (Fukazawa et al. 1994), showing evidence 
for such a gradient in the cooling region [r < rcooi ~ 6', or 60h~^kpc, assuming a pure Hubble flow 
with Hubble constant Hq = 100/ikms~-'^Mpc~''^). However, since the discrepancies between data 
taken with different instruments are larger than the formal error bars given, it is not clear what 
errors should be assigned to the data. Thus, we do not attempt a detailed fit, but only a qualitative 
comparison, and we also omit the error bars from our graphs. 

The Centaurus cluster has a bimodal redshift distribution (Lucey, Currie, & Dickens 1986), 
with two clumps, Cen 30 and Cen 45, each surrounding one of the two brightest elliptical galaxies, 
NGC 4696 {cz = 2923 km s'^) and NGC 4709 {cz = 4519 km s'^). The projected distance between 
these two galaxies is only 15'. The X-ray emission of the Centaurus cluster is clearly peaked at the 
position of the brightest galaxy, NGC 4696, with no evidence for a peak at NGC 4709 (Allen & 
Fabian 1994). 

Figure 1 shows the projected V-luminosity profile of Cen 30. NGC 4696 is approximated as a 
de Vaucouleurs profile with effective radius = 1.8' and total luminosity Ly = 8.1 x 10^'^h~^ Lq 
(Schombert 1987). The total V-luminosities of all other galaxies were taken from the data of 
Dickens, Currie, & Lucey (1986), using their conversions among wavelength bands and assuming 
de Vaucouleurs profiles for all galaxies. -"^ The luminosity distribution in and around the cooling flow 
region is strongly dominated by the central galaxy, even in projection, making it the likely source of 
any metallicity excess in this area. We take the metal injection function to be ([r) = r]i[r), where 
?7 is an adjustable parameter (metal injection per unit time per unit stellar luminosity), and i[r) is 
the V-luminosity per unit volume of the central galaxy, taken as a de Vaucouleurs profile with the 
parameters given above. 

Figure 2 shows the data from the ASCA Gas Imaging Spectrometer (CIS) and Solid-State 
Imaging Spectrometer (SIS; Fukazawa et al. 1994) and those from the ROSAT Position-Sensitive 
Proportional Counter (PSPC; Allen & Fabian 1994), together with results from different models. 
In all of the latter, we assume a power-law gas density profile with = 1.25, which gives a good 
approximation to the ROSAT X-ray surface brightness profile (Allen & Fabian 1994). The total mass 

^ This is not accurate in many cases, especially for spiral galaxies, but tends to overestimate the 
true luminosity of the galaxies. In any case, this correction is typically < 20%, and always < 41%. 
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distribution is not well known, since no lensing observations exist and the temperature profile is fairly 
uncertain (compare the results of Allen & Fabian 1994 with those of Fukazawaet al. 1994). Thus, we 
first consider two cooling-flow models which approximately span the range of observationally allowed 
temperature profiles. In the first, the mass profile, with Am = 1 + (3 — 2\p)/{l — a), is chosen to 
make M = constant, corresponding to a single-phase flow. In the second, we use an isothermal mass 
profile (Am = 1), which gives a multiphase flow with phase-dependent metallicity. Since our model 
assumes a steady-state cooling flow, which does not apply for r > rcooi, the metallicity differences 
among phases shown should be regarded more as indications of the size of the expected effect than 
as precise quantitative predictions. With the present measurement uncertainties, the disagreement 
among the data is somewhat larger than the expected iron abundance difference between different 
phases. This is partly because of the steep surface brightness profile of this cluster, which requires 
a rather narrow density distribution. 

In both models, the mass distribution is normalized by the emission-weighted temperature at 
the cooling radius determined (with low precision) by ASCA, Tg{rcooi) ~ 3keV. Above ~ 2keV, 
bremsstrahlung is the most important cooling mechanism, and the cooling function A(T) can be 
approximated by a power law. Below this temperature, A(T) increases dramatically due to line 
emission (e.g., Gehrels & Williams 1993). However, even if A(T) is approximated by a power law, 
the cooling time decreases substantially (oc T'^"") as the temperature decreases. Thus, the cooler 
phases cool quickly and nearly isobarically until they drop out of the flow, and their total emissivity 
is determined only by the mass deposition rate. Since the mass and volume of the cooling flow are 
dominated by the hotter phases, inferred properties such as M[r) will be unaffected by a change in 
the cooling function at lower temperatures, as long as the hotter phases remain above ~ 2keV.2 The 
metallicity of each phase at deposition is determined by its evolution at high temperatures rather 
than during the quick cooling through the low-temperature tail of the distribution, so it should be 
similarly unaffected by this change. Thus, a power-law cooling function should be appropriate for our 
purpose, namely, to calculate the metallicity as a function of radius, even though the high-density 
tail of the resulting density distribution is not correct. For these low temperatures, the cooling is 
practically isobaric and the emitted spectrum is independent of the structure of the cooling flow (see 
Johnstone et al. 1992). 

The shape of the metallicity profile is also quite insensitive to the value of a. For example, 
for M = constant, it is determined solely by C(^)i for the hottest phase in an isothermal mass 
profile, it depends on ([r) and A^. For definiteness, we choose a power-law cooling function with 
Ao = 7.8 X 10'^° erg cm^ g~'^ s~-'^ii'~°'^ and a = 0.2, which approximates the solar-metallicity curve 
of Gehrels & Williams (1993) near 3keV. 

In addition to these two standard cooling flow models, we consider two "toy" models with no 
mass deposition, i.e., in which a heat source is balancing the radiative cooling. In the first, we 
assume that the gas is static and accumulates locally injected iron over the lifetime of the cluster, so 
that Z[r) = Zoo + [•^(^)/p(^)] Jo v{t')d't'j where t is the age of the cluster in its present, relaxed state, 
and p is the local gas density. In the second, we assume perfect mixing within the cooling radius, 
due to a cycle of convection of the hot phases and evaporation of the cooled gas, as mentioned in 

This is not true at small r in the M = constant model, and therefore we expect this 
approximation to somewhat underestimate M and therefore the required 77 for this case. 
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§ 1. 

A visual inspection of the figures shows that the cooling flow models reproduce the observational 
curves somewhat better than the two extreme models with no net cooling. However, intermediate 
models with no cooling are possible, and even the extreme models may not yet be ruled out, 
given the uncertainties in the measured metallicities. The iron injection rates required by all the 
models considered (see the caption of Fig. 2) are in rough agreement with the observed rate of 
type la supernovae in elliptical galaxies (1 — 5 x 10~^^h^ yr""*^; Turatto, Cappellaro, & Benetti 
1994), assuming that each supernova injects a substantial fraction of a solar mass of iron (Nomoto, 
Thielemann, & Yokoi 1984). However, these rates are substantially higher than those inferred from 
x-ray observations of the hot gas around some elliptical galaxies not in cluster centers (Serlemitsos 
et al. 1993; Loewenstein et al. 1994). The expected injection rate from evolved stars is less certain, 
especially because of the uncertain calibration of stellar metallicity indicators (e.g., Worthey 1994). 
However, the metallicity of NGC 4696 is high compared to other elliptical galaxies (Davies et al. 
1987; Carollo, Danziger, & Buson 1993), implying that this process may well be significant. The 
relative importance of type la supernovae and stellar mass loss could be decided by spatially resolved 
abundance measurements of other metals. In particular, Fukazawa et al. (1994) have reported that 
O, Si, Ar, and Ca lines are also stronger in the central region, which would argue in favor of stellar 
mass loss. It should also be pointed out that a radial decrease in the stellar metallicity has been 
detected in the central part (~ 10") of NGC 4696 and other elliptical galaxies (Carollo et al. 1993). 
If this gradient continues out to arcminute scales, it could steepen the expected metallicity gradient 
in the gas. 



4. CONCLUSIONS 



The main conclusion we have reached in this paper is that the standard cooling flow model, 
where gas with the initial metallicity of the intracluster medium flows towards the center of the 
potential well, with the coolest phases continuously dropping out of the flow, is consistent with the 
metallicity gradient observed in the Centaurus cluster given reasonable rates of ejection of metals 
from Type la supernovae and evolved stars in the central cluster galaxy (NGC 4696). Most of the 
ejected metals are deposited in the hot phases, and these are the ones that survive for a longer time 
in the flow. The gas in the center of the flow is enriched owing to the metals that were deposited 
in the hot phases at larger radius. This model predicts that, at any radius, the metallicity should 
increase with the temperature of the phase. This predicts changes in the X-ray spectrum of cooling 
flows, which might be observable in the future. 

Alternative models would predict slightly different metallicity profiles. If the gas was static, 
the metallicity peak in the center should be more pronounced. In a model where a heating source 
produces a convection zone, where the gas follows a cycle where it is heated in the center, rises to 
large radius, and sinks again as it cools, the gas would be mixed and all the metals ejected would be 
present in the X-ray emitting gas, producing a central "plateau" in metallicity and a sharper decline 
at the edge of the convection zone. The metallicity would also be more uniform among the phases. 
However, in the absence of a specific model for a possible heating source, this alternative possibility 
is not sufficiently predictive to be tested from the metallicity profile. 



- 7 - 



The metallicity enrichment of the cooling flow gas depends mainly on the ratio of the stellar 
mass in the central galaxy, to the gas mass in the cooling flow. This ratio is relatively large in 
the Centaurus cluster, and the metallicity gradient should be much smaller in clusters with massive 
cooling flows where the central galaxy is not very luminous. At the same time, the metallicity 
gradient should depend on the merging history of clusters: it should be stronger in clusters which 
have grown quiescently by accreting many small clumps over a long period, without disturbing the 
radial distribution of the central gas, and weaker in clusters where the central gas has recently been 
stirred by a large merger. Thus, variations of the metallicity gradient among different clusters are 
expected. 
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FIGURE CAPTIONS 



Figure 1: Enclosed V-luminosity as a function of radius of the Cen 30 redshift clump (1700 km s""*^ < 
cz < 4150 km s"-*^) in the Centaurus cluster, based on the observations of Dickens et al. (1986) and 
Schombert (1987). The lines correspond to the luminosity of the central galaxy, NGC 4696 (solid), 
the early-type galaxies (dashed), and all galaxies (dot-dashed). 

Figure 2: Iron abundance as a function of radius in the region inside and surrounding NGC 4696. 
The points correspond to observations with three different instruments, the ROSAT PSPC (Allen 
& Fabian 1994), and the ASCA CIS and SIS (Fukazawa et al. 1994). The dashed curve is the 
prediction of the single-phase cooling flow model. The solid curves enclose the range of abundances 
in different phases in the multiphase cooling flow for an isothermal mass profile, M{r) oc r. The 
dot-dashed curve is the iron abundance accumulated during the age of the cluster by a static gas. 
For all three cases, the parameters Zoa and 77 = ?7_i3 x lO~-'^^M0i0''^ yr"-*^ were chosen to fit the 
SIS data, excluding the innermost point (which may be affected by ASCA's finite spatial resolution; 
Tanaka, Inoue, & Holt 1994). For the single-phase and M ex. r cooling flow models, ?7_i3 = 7.4 
and 3.2, respectively, and for the no-flow model, Hq J rj-isdt = 1.3. The dotted line corresponds to 
perfect mixing of the last model's iron mass within the cooling region. 
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